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In high-superconducting transition temperature ( T c ) iron-based superconductors, interband sign 
reversal (s±) and sign preserving (s++) s-wave superconducting states have been primarily discussed 
as the plausible superconducting mechanism. We study Co impurity scattering effects on the su¬ 
perconductivity in order to achieve an important clue on the pairing mechanism using single crystal 
Fei_ x Coa;Se and depict a phase diagram of a FeSe system. Both superconductivity and structural 
transition / orbital order are suppressed by the Co replacement on the Fe sites and disappear above 
x = 0.036. These correlated suppressions represent a common background physics behind these 
physical phenomena in the multiband Fermi surfaces of FeSe. By comparing experimental data and 
theories so far proposed, the suppression of T c against the residual resistivity is shown to be much 
weaker than that predicted in the case of a general sign reversal and a full gap s± models. The origin 
of the superconducting paring in FeSe is discussed in terms of its multiband electronic structure. 


I. INTRODUCTION 

The superconducting paring mechanism of high tem¬ 
perature superconductivity has been a long-lasting prior¬ 
ity research area, and is one of the most important and 
intriguing scientific subjects. After discovery of Fe-based 
superconductors (FeSCs), the superconducting mecha¬ 
nism has been discussed from the viewpoint of their 
unique multiband Fermi surfaces. The superconducting 
gap functions primarily discussed in FeSCs are the inter¬ 
band sign reversal s-wave (s±) and the sign preserving 
s-wave (s++) states. [TH3]. A stripe type antiferromag¬ 
netic spin fluctuation has been considered to mediate s±- 
wave state, while Fe 3 d orbital fluctuation mediates the 
s+_i_-wave state. 

Impurity scattering in superconductivity gives impor¬ 
tant information for understanding the pairing mecha¬ 
nism. A phonon mediated isotropic “BCS” superconduc¬ 
tivity is robust against the nonmagnetic impurity scat¬ 
tering. Meanwhile, since the Cooper pair is glued by the 
fc-dependent anisotropic scattering in the sign reversal su¬ 
perconductivity, e.g. ( 7 r, 7 r) spin fluctuation in cuprates, 
similar impurities induce random scattering by ending up 
with a strong pair breaking 015]. A general theory of 
pair breaking in the latter case was given by Abrikosov 
and Gorkov (AG-theory) 0. In the FeSCs, impurities 
are expected to induce a significant reduction in super¬ 
conducting transition temperature (T c ) for the s±-wave 
states. Although the nonmagnetic impurity doping ef¬ 
fects have been mostly examined by the AG-theory in 
connection to the scaling relation between T c and residual 
resistivity (po) jTUlOj ■ the important intrinsic multiband 
nature has been frequently neglected. 
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FeSe is one of the FeSC families showing superconduc¬ 
tivity at around 9 K mi. In the vicinity of the tetrag¬ 
onal / orthorhombic structural transition temperature 
(T s « 90 K), the orbital order has been recognized to 
develop without long-range antiferromagnetic (AFM) or¬ 
der, being in strong contrast with the other FeSCs [12b 
£7]. A sign preserved superconducting state has been 
indicated by scanning tunneling microscopy (STM) and 
spectroscopy (STS) in single layer FeSe [H]. Neverthe¬ 
less, in thicker films and bulk FeSe, the nodal supercon¬ 
ducting gap has been suggested by STM/STS as well 
as the London penetration depth, implying the contri¬ 
bution of AFM spin fluctuation to the formation of the 
Cooper pairs HSIED]. Both neutron scattering and nu¬ 
clear magnetic resonance have shown AFM spin fluc¬ 
tuations mum [mug, whereas an imperfect nesting 
between electron and hole Fermi surfaces has been ob¬ 
served by orbital resolved angle resolved photo emission 
spectroscopy (ARPES) [23] ■ Detailed studies on Co dop¬ 
ing are considered to give an important hint for under¬ 
standing the mechanism of superconductivity presently 
debated in FeSe. 

In the present paper, we report systematic electrical 
transport studies on the effect of Co impurity doping for 
single crystal FeSe. In the FeAs system, Co has been 
regarded as a nonmagnetic impurity and provides an ad¬ 
ditional electron as a carrier [21], being in contrast with 
the situation of Mn ,25], [26]. Present Hall coefficients 
and ARPES [27] [2B] support this understanding. There¬ 
fore, Co is considered to act as a nonmagnetic impurity 
and add an electron as an itinerant carrier to FeSe. Our 
present experimental data indicate a correlation between 
the suppression of T c and that of the structural transition 
orbital order when Fe is replaced by Co. By carefully an¬ 
alyzing the dependencies between the suppression of T c 
and po and taking into consideration the realistic elec¬ 
tron and hole Fermi surfaces in FeSe [22], we found that 
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FIG. 1. (a) The X-ray diffraction data of Fei-sCo^Se (i n0 m = 0, 0.010, 0.025, 0.05, and 0.10) around the c-axis. The inset 

represents the normalized (004) peaks, (b) The nominal concentration (x n om) dependence determined by EDS (zeds) and 
Tkink- Tkink is defined as the temperature showing a kink in the resistivity curve as described in the later section, (c) The 
nominal concentration dependence of concentration Fe+Co with respect to Se (i.e. (Fe+Co)/Se). The error bars are defined 
as the standard deviations of measurements. 


the suppression of T c against po is clearly weaker than 
those expected from both the AG-theory and a recent 
more particular model for the s±-wave state [30]. Our ex¬ 
perimental observations give better agreement with the 
sign preserved superconducting gap states and suggest 
important multiorbital nature of Fermi surfaces for the 
occurrence in superconductivity of FeSe [31 GEM- 


II. EXPERIMENTS 


High quality single crystals of Fei-^Co^Se (0 < 
x < 0.075) were grown by a molten salt flux method 
pH [331 [34]. The temperature of hot and cold positions 
was kept for 390 and 240°C, respectively. After ~ 10 
days, single crystals were grown around the cold part of 
the quartz tube. Being different from the conventional 
method, polycrystalline samples |35j were employed as a 
precursor. Note that z nom indicates the nominal compo¬ 
sition of Co applied for synthesis of polycrystals. We 
examined the quality of the prepared samples by X- 
ray diffraction (XRD, Cu K -a radiation wavelength of 
1.5406 A) around the c-axis as well as energy dispersive 
X-ray spectroscopy (EDS). The magnetic susceptibility 
was measured at B = IT parallel to the ab-plane. The 
temperature dependence of electrical resistivity and Hall 
resistance was measured by a four probe method. The su¬ 
perconducting transition temperatures (T c s) were deter¬ 
mined at the end point of the superconducting transition 
with a value of approximately less than l.OxlO -8 Hem. 
In the Hall resistance measurements, magnetic fields were 
varied in a range of \B\ < 9 T parallel to the c-axis. 


III. RESULTS 
A. Sample Characterization 

Figure [I] (a) shows XRD spectra obtained for i n om 
= 0, 0.01, 0.025, 0.05, and 0.10. Clear (00Z) (Z = 1-4) 
peaks were observed. No impurity peaks were detected 
within the experimental errors. Since the lattice shrink¬ 
age is very small in the case of Co doping [35], it was 
difficult to see a significant influence on their lattice pa¬ 
rameters among samples with small concentration of Co 
even though the shift in the (004) peak between z nom 
= 0 and 0.10 was detected. The black squares in Fig. 
[l] (b) represent the dependence of z nom as a function 
of Co concentration determined by EDS (zeds)- EDS 
spectra were taken at 10 different positions on each spec¬ 
imen. The error bars were estimated as the standard 
deviations from the average. A monotonic increase in 
zeds was observed against z nom . Note that the nominal 
concentration of z nom = 0.01 is used for the sample in¬ 
stead of Zeds because EDS was not sufficiently sensitive 
for detecting the Co peak under such a small concentra¬ 
tion. The blue circles in Fig. |T] (b) show that the Tkink 
(the kink temperature in the resistivity curve as shown 
in the section of the electrical resistivity) reduced mono- 
tonically with z nom . These experimental data indicate 
that the systematic substitution of Co was successfully 
accomplished. We regard zeds as the real x and used 
this for discussion hereafter. In the parent FeSe grown 
by a KCI/AICI 3 flux method, no interstitial Fe was found 
by X-ray diffraction measurements |33j . In order to check 
the interstitial Fe carefully in Fei-^Co^Se, the concentra¬ 
tion ratio: (Fe+Co)/Se was compared against the nomi¬ 
nal Co concentrations as shown in Fig. [l] (c). Since the 
ratio is almost constant within the error bars for all nomi¬ 
nal concentrations, the interstitial Fe is not influenced by 
the Co doping. It should be emphasized that we selected 
the samples with (Fe+Co)/Se < 1.05 for electrical trans- 
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port measurements. 

B. Electrical Resistivity 



FIG. 2. (a) Normalized temperature dependence of resistivity 
(p/p300K). Note that each curve is shifted by 0.1 for clarity. 
The closed triangles indicate the kink in the resistivity curve 
(Tkink). Tkink is defined as the peak position of the 1st T 
derivative of p. (b) Temperature dependence of the 1st deriva¬ 
tive of p. (c,d) Enlarged view of temperature dependence of 
p. The curve is not shifted, (c) The black and the red lines 
show fitting results of p-T curves in the normal states using p 
= (aT+b) or p = (aT 7 +/3), respectively, (d) Closed and open 
black circles represent x = 0. The others (blue, yellow, and 
green) indicate x nom = 0.025. Their averaged concentration 
was x ~ 0.018. 

Temperature dependence of the normalized electrical 
resistivity (p/p 3 ook) for Fei-^Co^Se (0 < x < 0.075) is 
shown in Fig. [2ja). In parent FeSe, the p/p 3 ook showed a 
kink (T kink ) due to the structural transition [TT}[33|. This 
can clearly be seen as a dip of the first derivative as shown 
in Fig. [2])b). T k j nk decreased with an increase in Co 
concentration and finally disappeared at x = 0.036, indi¬ 
cating that the high-temperature tetragonal phase holds 
above x = 0.036. As seen in the magnified p-T plots (Fig. 
[2jc) ) at low temperatures, the T c decreases with an in¬ 
crease in Co and the superconductivity is killed at x = 
0.036 above 2 K. Both suppression of superconductivity 
and structural transition in FeSe occur by partial replace¬ 
ment of Fe to Co. In order to interpret the T c suppression 
by Co in terms of the impurity scattering, residual resis¬ 
tivity po was evaluated. The absolute values of resistiv¬ 
ity (p) as a function of temperature are very sensitive to 


the measurement conditions and generally contain large 
errors arising from sample shape or/and inhomogeneity, 
etc.. In order to estimate a reliable value, resistivity mea¬ 
surements were carefully performed on more than twenty 
samples selected from the same sample batch of * n0 m = 
0.025. Figure [2](d) shows the results by excluding the 
samples showing large deviation from the average. Just 
above the T c s, p of a; nom = 0.025 (x « 0.018) are clearly 
larger than those of FeSe. 

In order to deduce residual resistivity (po), two types 
of fitting were employed. One is linear fitting using a 
function of aT + b , where b represents po- The other 
is fitting using a power function of aT 7 + /3, where (3 
represents po- The fitting was made in the temperature 
range of 12^13K (just above the superconducting onsets) 
to 25~27 K. An example of fitting line is shown in Fig. 
i< 0 - Figure [3] shows the x dependence of these fitting 
parameters. Although some errors can be recognized, po 
apparently increases monotonically with x as seen in Fig. 
[3])a). The value of po was not significantly affected by the 
fitting functions. The exponent 7 in the power function is 
depicted against x in Fig. [3] (b). At small concentrations 
it slightly decreases from ~ 1.1 to ^ 0.8 with x, and then 
jumped to 7 « 2.0 at x ~ 0.08. 
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FIG. 3. (a) Co concentration ( x ) dependence of the residual 
resistivity, i.e. b in the aT + b and /3 in the aT 7 + /?, versus 
x. (b) The exponent 7 versus x. 


C. Magnetic Susceptibility 

Temperature dependence of magnetic susceptibility 
( M/H ) was measured on field cooling for Fei-^Co^Se 
(x ~ 0, 0.018, 0.034) as shown in Fig. [4j A kink was 
observed in the temperature dependence of M/H curve 
at the structural transition temperature in FeSe, being 
in consistent with the earlier study [33]. We express this 
temperature as T mag in the present paper. The T mag , 
the transition temperature in spin magnetic moments, 
was also found in our samples of x ss 0 , 0.018, but it 
disappears in the sample of x « 0.034. This disappear¬ 
ance, most presumably corresponding to the absence of 
the tetragonal to orthorhombic structural transition, co¬ 
incides with the disappearance of T k i n k and T*. 

In order to confirm the mangetic statement of Co, 
M/H curves were fitted in a low temperature range (15 
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FIG. 4. Temperature dependence of magnetic susceptibility 
( M/H ) for Fei-zCo^Se with (a-c) * R5 0, 0.018, 0.034. (b,c) 
Enlarged views are inserted. B = IT was applied in paral¬ 
lel to the afe-plane of crystals. All data were taken on field 
cooling. The Tmag is defined as the kink of the curves. The 
red lines indicate Curie-Weiss fitting at low temperature, i.e. 
C/T + A where C and A are the Curie constant and a certain 
real valued constant, respectively. The black lines also repre¬ 
sent the Curie-Weiss curves with C s i m simulated by assuming 
that Co is a magnetic impurity with the angular momentum 
quantum number J = 1/2. 


-60 K) by a function C/T + A , where the Curie con¬ 
stant is assumed to be C and A is a certain constant 
value. As shown in Fig. [4j the fitting curves are repre¬ 
sented as red lines and the obtained C merely changes 
as a small value with an increase in x. For comparison, 
the Curie constants were simulated (Csim) by assuming 
that all Co elements are localized with spin 1/2. The 
evaluated black curves show large deviations from the ex¬ 
perimental plots. The obtained values of Curie constant 
were much smaller than the doped cobalt concentration 
and did not change significantly. This result supports 
the scenario that the cobalt is acting as a nonmagnetic 
impurity in the Fei_a,Co x Se system. 


D. Hall Resistivity 

The magnetic field ( B ) dependence of Hall resistivity 
( p yx ) is shown in Fig. [5] In order to remove the extrin¬ 
sic deviations from symmetric components caused by the 
misalignment of electrodes, corrections were made on the 


raw p yx curves: 


m + i in 

A strong nonlinear behavior at low temperatures was 
moderated as x increases. This can be understood to 
be caused by the suppression of mobility. Moreover, 
temperature dependence was also moderated, being con¬ 
sistent with Hall coefficient’s (Ah’s) T dependence de¬ 
scribed later. Ah was determined from the slope of p yx 
at low B where p yx s develop linearly. The linear fitting 
is displayed as a dashed line in Fig. [11(a)- 

Temperature dependence of Ah is shown in Fig. [6] 
For the parent FeSe, Ah became positive at around Tki„k 
and dropped with decreasing temperature. We define T* 
as the temperature where Ah shows a sudden decrease. 
For x = 0.017, T* is derived from the intersection of two 
linear fitting lines as shown in the gray lines of Fig. [6] 
With an increase in Co concentration, T* lowered and 
almost disappeared at x = 0.036, being consistent with 
the disappearance of Tkmk- The amplitude of Ah at low 
temperatures once increased and then decreased with an 
increase in x. When the concentration of Co exceeded 
above x = 0.036, the observed change in Ah is nearly 
suppressed. 


IV. ANALYSES OF T c SUPPRESSION RATE BY 
CO DOPING 

Figure [7] shows the x dependence of T c . T c monoton- 
ically decreased with an increase in x within the exper¬ 
imental errors. Figure [8] represents the T c as a function 
of po . The po values were estimated from the intercepts 
of the linear fitting (white diamonds) and the T 7 fitting 
(blue circles) on the resistivity curves at low tempera¬ 
tures as shown in Fig. [2](c). T c monotonically decreased 
with an increase in po and disappeared at po of « 60 - 80 
pflcm. 


A. AG-theory in realistic two band model 

In order to clarify the superconducting pairing mech¬ 
anism in FeSe, here we evaluate how T c is suppressed 
against residual resistivity when Co is doped as an im¬ 
purity in the light of two models. The first model is 
the AG-theory which can be applicable to the sign rever¬ 
sal superconducting states [6]. The other is the specific 
model formulated for the full gap s±-wave states theoret¬ 
ically calculated based on the realistic physical condition 
of FeSCs [30] . Both models are valid in the limit that 
nonmagnetic impurities are not incorporated to a great 
extent. The AG-theory is formulated as, 
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FIG. 5. Magnetic field ( B ) dependence of Hall resistivity ( p yx ) of x = 0, 0.010, 0.017, 0.036, and 0.075 ((a)-(e), respectively). 
The linear fitting line to derive the Hall coefficient (.Rh) in the low B limit is shown as a dashed line in (a). The p yx s were 
measured at different temperatures for x = 0 as noted in the legend at bottom left. 



T( K) 

FIG. 6. Temperature dependence of Hall coefficient (Ah) es¬ 
timated in the low magnetic field regime. The black arrows 
indite the temperature where Ah starts suddenly decreasing 
(T*). For x = 0.017, T* is defined as the intersection point 
of two linear fitting lines shown as gray in color. 


where /cb is the Boltzmann constant, T c ° is the supercon¬ 
ducting transition temperature in the pristine material, r 
is the scattering time, and 'P(z) is the digamma function. 
Here we use superconducting transition temperature as 
T^ G to differentiate from experimental values ( T c ). The 
pair breaking energy is characterized by r. The value 
of r was computed at each value of T^ G using Eq. [ 2 ] 
with T c0 = 9 K. In order to compare the AG-theory and 
our experimental data, r needs to be converted into the 
residual resistivity in a reasonable way as follows. 

In FeSe, nearly semimetallic electronic structure is re¬ 
alized [Til rra I2U1IMI55] . The ARPES reported a single 
electron Fermi surface (FS) around the M-point and a 
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FIG. 7. The Co concentration ( x ) dependence of the super¬ 
conducting transition temperature ( T c ). Down arrows indi¬ 
cate that the zero resistivity is not observed above 2 K. The 
gray broken line is guide to the eyes. 

single hole FS around T. It should be mentioned that 
the recent mobility spectrum analysis and the quantum 
oscillation measurements in magnetoresistance reported 
an additional small FS El EZ]. In the present work, 
since an influence of this tiny FS may be small in the 
present analysis, we focus on the two main FSs observed 
by ARPES. In this system, therefore, conductivity can 
be described in a nearly free electron model by keeping 
the effective mass (m*) constant on each FS as follows: 

cr = 7 7l e 2 n+ry 2 e 2 r 2 . (3) 

where r)i = rii/m*{i=1 ,2) and rii is a carrier density of 
*-th FS. From this equation we define a weighted average 
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T (rave) as, 


Therefore, 


7~ave 


Vl Tj + ?72T~2 
??1 +??2 


(T = +77 2 )e 2 r a ve- 


(4) 

(5) 


Thus, the relation between the residual resistivity (p AG = 
1/a) and T ave , which should correspond to the r in the 
Eq[2] is derived. The information on FeSe achieved from 
the ARPES measured at 30 K to calculate r)i is shown in 
table I pTj [29]. Since the anisotropy of the hole FS at T 
is smaller than the electron FS at M-point, the observed 
ARPES may represent such an averaged value. For each 
FS, m* was calculated by the relation m* = H 2 k f/v^ 6 . 
The vf was averaged over the electron FS and was eval¬ 
uated as i>p ve . The general carrier density definition is, 


n = 


2 

J2n? 



( 6 ) 


where the factor of 2 is for the spin degeneracy and the 
integration is performed for the entire volume surrounded 
by FS. Given the elliptic cylinder, n can be calculated as, 


n = 


U ;>s 

/bp/bp 

27TC 


(7) 


where kp and kp denote the Fermi wave number of the 
longer and the shorter axes of the ellipse, respectively. 
The c-axis cell parameter expressed as c in the real 
space is taken from literature [38]. Now we can calcu¬ 
late the p AG in the framework of nearly free electron 
approximation from Eq. [5] The obtained parameters 
are ni ss 2.90 x 10 26 nr - , ?z 2 ~ 3.35 x 10 26 m -3 and 
m\ ss 2.28 x 10 -30 kg, m/, ~ 4.26 x 10 _3 ° kg, where the 
indices of 1 and 2 represent the electron and the hole 
bands, respectively. 

The obtained dependence, i.e. T C AG versus p AG , is rep¬ 
resented by a solid red line in Fig. [8] The suppression of 
T ag is clearly faster than our experimental data. Given 
a conventional single-band model, for comparison, the 
analytical expectation (dashed line) is represented in the 
same figure, where the carrier number n was evaluated 
from i?H data measured at low B’s. The averaged ef¬ 
fective mass was employed to calculate p AG . The sup¬ 
pression of T AG is faster in the single-band model than 
that in the two-band model. Since the amplitude of J?h 
is suppressed in a compensated semimetallic state, the 
single-band model (i?H = 1/ne) overestimates n and con¬ 
sequently leads to faster suppression of T AG . 


Electron doping effects on the two band model 

We interpret the electron doping effects on our fitting 
model based on a nearly free electron model. If the carrier 


TABLE I. Results of the ARPES measurements at 30 K [14 , 
l29j . Note that several photon energy results are averaged. 


r 

M (long axis) 

M (short axis) 

fc F (A -1 ) 0.11 
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FIG. 8. Dependence of the superconducting transition tem¬ 
perature ( T c ) on the residual resistivity (po). The po was esti¬ 
mated by fitting p-T curves in the normal states as shown in 
Fig.® c). White diamonds and blue circles represent po esti¬ 
mated by T-linear and T 7 type fitting functions, respectively. 
With realistic band parameters, the red, blue and cyan lines 
are drawn by the AG-theory (T C AG versus p AG ) and the s±- 
wave model (T c s± versus p^ ) with and without orbital order 
(00), respectively. The dashed line represents the AG-theory 
in the conventional single carrier model. 


doping rigidly shifts the Fermi energy (Ef) as represented 
in Fig. [9] the conductivity is shown not to be affected 
by Vi + V 2 but determined only by r ave as follows: In a 
nearly free electron model, energy dispersion relation of 
an electron band is written as, 

^ < 8 > 

The amplitude of Fermi wave number is derived by 
putting E{ k) = 0 as |kp| = dc\/2m* Ep /H. The carrier 
density of a 2 dimensional FS (a cylinder) is, 


|kp| 2 m*Ep 
27rc nh 2 c 


(9) 


By assuming electron doping by 5, Ef becomes Ef + 6. 
The value of in* is not modified in a rigid band model. 
Therefore the conductivity after the electron doping can 
be written as, 
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a = (77! + r ]' 2 ) e 2 T' ve 

= ( eP+6 E™-6 \ , 

^ 7Tfi 2 C 7 tS 2 C J ^ 

= (771 + m )e 2 T' ve (10) 

where .Bp 1 -* and Ep 2 ^ represents the Fermi energy of elec¬ 
tron and hole bands, respectively. Therefore, (771 + 772 ) 
of a 2 dimensional nearly free electron model does not 
change by rigid band shift. As long as the approximation 
is valid, the effect of Co substitution can be simplified to 
the scattering time problem. 


(a) 




FIG. 9. A schematic view of Fermi surfaces of (a) two band 
semimetallic structure and (b) electron doped one on the k x — 
k y plane. The bottom figure represents the rigid band Fermi 
energy shift due to the electron doping by 8. 


B. Theory of impurity-induced reduction in T c for 
the full gap s±-wave state 

The second model to analyze the suppression of T c is a 
realistic multiorbital theory. Here, we employed a multi¬ 
orbital tight-binding model for FeSe with realistic Fermi 
surfaces in order to perform quantitative theoretical stud¬ 
ies of the impurity effect on the s±-wave state. Follow¬ 
ing two theoretical models are introduced: In model (a), 
the orbital order is absent and the C\ symmetry is pre¬ 
served. In model (b), the experimentally observed sign- 
reversing orbital polarization is introduced, by which the 
symmetry is lowered to C 2 . The Fermi surfaces of these 
models without carrier doping are shown in Figs. [TO] 
(a) and (b), respectively. The carrier numbers per Fe 
site are rq, = 2 . 02 % and n e = 2 . 02 % in model (a), and 
rih = 1.95% and n e = 1.94% in model (b). In model 
(a) (model (b)), the hole-pockets disappear by introduc¬ 
ing 3.0% (3.7%) electron carriers induced by Co-doping. 
This result would be consistent with the disappearance of 
the structure transition by 3.6% Co-doping as observed 
experimentally. 

In order to construct the model (a) without the or¬ 
bital order, we first performed LDA band calculations 


for FeSe using the wien2k package, and made the tight- 
binding fitting using the wien2wannier package. How¬ 
ever, the obtained Fermi pockets are too large. To repro¬ 
duce the tiny Fermi pockets experimentally observed in 
FeSe, we introduced the orbital-dependent energy shifts 
at r, X(Y), and M points; (5Ey,8E x ,8Em)'- We set 
( 8E t ,8E x ,y ,8E M ) = (—0.26,0.13,0) for xz,yz- orbitals, 
and (8Et,8Ex,y,8Em) = (0,0.26,-0.25) for cry-orbitals 
in the unit of eV. Such level shifts were introduced by 
the additional intra-orbital hopping integrals; 8t on ~ slte = 
SE r / 4 + SE x ,y/ 2 + <5£ m /4, 8t nn = 8E r / 8 - 8E x ,y/ 8 , 
and St nnn = SEr/16 + 8 E x ,y / 16 — 5Em/8 . In the model 
(b), we further introduced orbital polarization with sign- 
reversal: E yz — E xz = —0.02 eV at T and M points, and 
E yz — E xz = 0.02 eV at X and Y points. 

Figure 10 (c) shows the impurity-induced reduction in 
T c of the s±-wave state for T c 0 = 9 K as a function 
of residual resistivity po , obtained by applying the T- 
matrix approximation. We used the realistic first prin¬ 
ciples impurity potential model for Co-ion given in !39j . 
and neglected the carrier doping effect by Co doping. 
Here, we introduced the effective intra-orbital and inter¬ 
pocket repulsive pairing interactions due to the spin fluc¬ 
tuations, by adjusting the strength to realize T c0 = 9 
K. The detailed explanations for the gap equation anal¬ 
ysis had been given in |30j . In the spin-fluctuation pair¬ 
ing mechanism in FeSe, only the electrons with xz and 
yz orbital characters contribute to the superconductiv¬ 
ity since the a;y-orbital is absent in the hole-pocket HD] , 
We also put the renormalization for the xz, yz-orbitals 
(z) to be unity. In both models (a) and (b), the super¬ 
conducting state disappears when the residual resistivity 
is Pq 1 ss 4[yHcm], insensitively to the orbital polariza¬ 
tion. The current vertex correction for po was taken into 
account, and po was kept independent of z. If the ex¬ 
perimental value z ~ 1/3 for the xz, yz-orbitals is taken 
into account m, the critical residual resistivity is mul¬ 
tiplied by z~ x . Therefore, p$ m z _1 4[/if2cm]. Note that 
the total carrier number 7 t carr i er is ~4% per Fe ion in 
both model (a) and model (b). In a model with 77 C arrier 
~2.5% by reducing the sizes of Fermi pockets, the rela¬ 
tion pg r « 2 “ 1 6 . 5 [/xHcm] is realized. 

Based on the p$ ~ z _ 1 4[/rflcm] with z _1 ~3, the su¬ 
perconducting transition temperature (T c s± ) suppression 
predicted in the s ± -wave state is displayed by the blue 
(the orbital order is considered) and cyan (the orbital or¬ 
der is not considered) lines in Fig. [ 8 ] In both cases, the 
calculated results show faster reductions of T c s± than the 
experimental data, being consistent with that in the first 
model. 


V. DISCUSSION 

A. Electronic Phase diagram 

The experimental observations in Fei-^Co^Se are 
summarized in a phase diagram shown in Fig. H3 The 
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FIG. 11. An electronic phase diagram of Fei-^Co^Se. The 
superconducting transition temperature (T c : red circles), the 
temperature showing a kink in the p-T curve (Tkmk: green 
squares), and temperature showing a sudden decrease in the 
Hall coefficient (T*: light blue triangles) are represented. 
T mag (blue diamonds) represents a kink observed in temper¬ 
ature dependence of magnetic susceptibility. 


FIG. 10. (color online) Fermi surfaces of the present tight- 
binding models for FeSe. (a) without orbital polarization and 
(b) with orbital polarization with sign reversal between T and 
X(Y) points, (c) Obtained impurity-induced reduction in T c 
for the s±-wave state as a function of residual resistivity po 
(T c = 9 K), in both models (a) and (b). 

i?H at low temperatures increases first and then decreases 
with an increase in x of Co as shown in Fig. [6] Above 
x = 0.036, i?H becomes nearly constant against x. In a 
two band semi-classical model, the low field limit of Hall 
coefficient can be written as, 

_ pv 2 -np 2 
Ah , . \9 > 

eypv + npL)~ 

where n, p and /r, v are the carrier concentrations and the 
mobilities of holes and electrons, respectively. ARPES, 
STM and transport studies have demonstrated that the 
electronic structure in the orthorhombic FeSe is almost 
compensated semimetal m mi m m- In a compen¬ 
sated semimetal (n ~ p), the amplitude of Ah is sup¬ 
pressed due to the electron and hole compensation as 
shown in eq. [TT] When electrons are doped to FeSe, the 
compensation breaks to increase the amplitude of Ah- 
Further electron doping may shrink and vanish the tiny 
hole pocket. In such a case, since a single carrier picture 
is more appropriate, the amplitude of Ah is suppressed. 
Therefore, our experimental observations on the Co dop¬ 
ing dependence Ah in Fei-xCo^Se can qualitatively be 


understood in terms of the electron doping by Co. For x 
= 0.036 and 0.075, where the high temperature tetrago¬ 
nal phase remains in an entire temperature range, the Ah 
at 10 K gives the carrier number (n) of ~ 10 21 cm -3 in 
a single carrier model (i.e. Ah = 1/n). By assuming the 
rigid band shift of the Fermi energy in Fei_ x Co x Se, the 
carrier number in the electron pocket above x = 0.036 is 
deduced to be ~ 10 21 cm -3 based on the carrier numbers 
of the electron and the hole pockets in FeSe 23) being 
comparable with those in x = 0.036 and 0.075. Since the 
variation observed in Ah is consistent with the electron 
doping provided by Co, potential scattering and changes 
in both the type and the number of carriers influence on 
the phase diagram in Fig. [TTJ Note that quantitative 
analysis centering on the Hall resistivity curvature is not 
realistic in the present magnetic field regime. The trans¬ 
verse magnetoresistance is ss 10% at T =10 K and B = 9 
T for x ~ 0.018 whereas it is ss 300% in the parent FeSe 
[.Til , suggesting the significant suppression in carrier mo¬ 
bility by the Co scattering effects. Since only a limited 
number of carriers contribute to the magnetotransport 
below 9 T, any reasonable analysis on the magnetotrans¬ 
port data would not accurately be valid [J5] H3] . Semi- 
classical multicarrier fitting of Hall resistivity assuming 
the number of electron and hole pockets could not give 
correct transport parameters for anisotropic Fermi sur¬ 
faces, and somewhat more sophisticated analyses such 
as a mobility spectrum method are necessary [34] . For 
preventing these uncertain factors in discussion, we em¬ 
ployed the band parameters extracted from the ARPES 
data for the present discussion. 














9 


Both superconductivity and structural transition are 
suppressed by Co doping and disappear at x = 0.036. 
This is consistent with the disappearance of the kink in 
magnetic susceptibility curves (T mag ). Since the orbital 
order was reported to develop in the orthorhombic phase 
of FeSe, the present results may also indicate the suppres¬ 
sion to be correlated between superconductivity and or¬ 
bital order by Co doping. Theoretically, nesting between 
electron and hole pockets is predicted to be important for 
both superconductivity and orbital order [l]-j3] in FeSCs 
32, 44, 45. In FeSe, electron doping may also suppress 
the nesting due to the small Fermi energies [2D], being 
consistent with the experimentally observed evolution of 
the Ru as a function of x. Hence, the suppression corre¬ 
lated among superconductivity, structural phase transi¬ 
tion, and orbital order may originate from the reduction 
in nesting between electron and hole pockets. 

B. Paring Mechanism 

If the sign reversal superconducting mechanism could 
be dominant in FeSe, suppression in superconductivity 
would be similar to that predicted by both theories of AG 
and the full gap s±-wave states. However, the suppres¬ 
sions experimentally observed as T c was much smaller 
than those described in these models. Consequently, the 
sign preserved superconducting mechanism will be more 
appropriate for Fei-^Co^Se. In FeSe, a nodal super¬ 
conducting gap was experimentally observed in the elec¬ 
tron pocket EDI®. Meanwhile, orbital resolved ARPES 
demonstrated an imperfect antiferromagnetic nesting be¬ 
tween the electron and the hole Fermi surfaces. One can 
imagine that the sign reversal and the sign preserved pair¬ 
ing mechanism may compete with each other via the in¬ 
ter and the intra orbital scattering HH31EU US- In this 
case, the sign preserved pairing mechanism may give a 
primary contribution in FeSe from the viewpoint of the 
T c suppression observed in our experimental data. An¬ 
other possible scenario may be that the superconducting 
gap structure changes from a nodal s- to an s++-wave 
states. It has been pointed out that the node of the su¬ 
perconducting gap is not protected by symmetry, and an 
accidental node ESI S3 was suggested. The supercon¬ 


ducting gap may change from a nodal gap to a full gap 
state depending on the quality of single crystals. In such 
a case, the potential scattering induced by Co may lift 
the node to open a superconducting gap, resulting in a 
change of gap structure from the nodal- to the s++-wave 
states. It should be noted that sign preserved super¬ 
conducting gap structure has been proposed for single¬ 
layered FeSe on SrTiC >3 [TB] possessing highly electron 
doped band structure [JSj. Our observations are consis¬ 
tent with these previous works. 

VI. CONCLUSION 

We investigated the nonmagnetic Co impurity doping 
effect on the superconductivity and the structural tran¬ 
sition/orbital order in FeSe. In the electronic phase dia¬ 
gram of Fei-^Co^Se, both superconductivity and struc¬ 
tural transition/orbital order were suppressed by Co dop¬ 
ing and disappeared at x = 0.036. Experimental data of 
the T c suppression against the residual resistivity were 
compared to the values expected from the AG theory [6] 
and a full gap s± theory [30]. Our data showed much 
slower suppression than those from both theories. As a 
result, the sign preserved superconducting mechanism is 
more appropriate to explain the superconducting pairing 
mechanism in FeSe. In FeSe, the interband Fermi surface 
nesting including both intra- and inter-orbital scattering 
processes play an essential role for the superconducting 
pairing mechanism and for the normal state phase dia¬ 
gram [32 IHJ 05] . The present studies unveiled a part of 
the complex electronic states. Considering the electronic 
structures taking place in this system, small perturba¬ 
tion may change the situation more complex in the Fe-3d 
multiband system by changing one state to another. 
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